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FINAL  REPORT 


A  Research  Program  on  Ultrawideband  Electromagnetic  Pulse 
Propagation  Through  Causal  Media 

Dr.  Natalie  A.  Cartwright 
State  University  of  New  York  at  New  Paltz 


The  general  purpose  of  AFOSR  grant  “A  Research  Program  on  Ultrawideband  Elec¬ 
tromagnetic  Pulse  Propagation  Through  Causal  Media”  (FA9550-13-1-0013)  was  to 
obtain  a  deeper  understanding  of  electromagnetic  pulse  propagation  through  disper¬ 
sive  material  and  subsequently  use  this  understanding  to  look  for  new  ways  to  improve 
technology.  The  significant  work  accomplished  through  this  grant  was  in  the  area  of 
waveform  design  for  synthetic  aperture  radar  imaging  through  dispersive  material. 


Introduction.  When  an  electromagnetic  pulse  travels  through  a  dispersive  material 
each  frequency  of  the  transmitted  pulse  changes  in  both  amplitude  and  phase,  and 
each  frequency  at  its  own  rate.  As  a  consequence,  broadband  pulses  propagating  in 
dispersive  material  experience  significant  amplitude  distortion  and  changes  in  pulse 
velocity.  Asymptotic  analysis  of  the  exact  integral  representation  of  the  propagated 
field,  which  utilizes  the  full  causal  dispersion  relation  of  the  dispersive  material,  pro¬ 
vides  a  complete  far-field  description  of  the  propagated  pulse.  In  1975,  Oughstun  and 
Sherman  [1]  utilized  asymptotic  methods  to  show  that  in  a  dielectric  (nonconducting) 
material  the  low-frequency  component  of  the  propagated  field,  the  so-called  Brillouin 
precursor,  has  a  peak  amplitude  that  decays  algebraically  with  propagation  distance 
(as  the  inverse  square  root  of  propagation  distance) ,  whereas  other  pulse  components 
decay  exponentially.  In  2005,  Oughstun  [2]  proposed  a  Brillouin  pulse  as  the  pulse 
that  experiences  the  least  attenuation  due  to  material  absorption.  It  has  been  hy¬ 
pothesized  that  the  slow  decay  rate  of  the  Brillouin  pulse  can  be  used  to  advantage 
in  radar  detection  and  imaging  applications  [2,  3,  4]. 

The  question  of  whether  or  not  the  Brillouin  precursor  is  the  ideal  waveform  for 
synthetic  aperture  radar  imaging  through  dispersive  material  was  partially  addressed 
by  Varslot,  Morales,  and  Cheney  [5,  6]  in  a  pair  of  papers  appearing  in  2010  and  2011. 
These  authors  used  a  filtered  back-projection  algorithm  to  derive  an  optimal  filter 
and  its  associated  optimal  waveform.  The  optimization  is  based  upon  minimizing 
the  mean  square  error  of  the  L2-norm  between  the  ideal  image  and  the  reconstructed 
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image.  The  optimal  waveform  is  numerically  derived  for  each  noise  level  by  solving 
a  fixed  point  equation  for  the  absolute  value  of  the  optimal  pulse  spectrum  and  then 
using  a  minimum-phase  algorithm  to  obtain  the  optimal  waveform.  The  authors 
concluded  that  the  optimal  minimum-phase  waveforms  have  a  “transmit  spectrum 
that  is  concentrated  around  the  frequencies  which  are  conducive  to  the  generation  of 
precursors”  [6],  but  no  conclusive  statement  could  be  made. 

The  results  of  our  research  effort  provide  an  asymptotic  expansion  of  the  propagated 
electric  field  component  of  the  impulse  response  due  to  scattering  by  an  isotropic 
point  source  in  a  frequency-dependent  dispersive  (and  lossy)  material.  We  define  a 
scattering  precursor  based  on  this  asymptotic  expansion.  The  scattering  precursor  is 
analogous  to  Oughstun’s  Brillouin  pulse  [2]  in  that  it  experiences  near  optimal  (if  not 
optimal)  penetration  into  the  material. 


Formulation.  We  adopt  the  assumptions  and  formulation  of  [5,  6],  one  of  which 
assumes  the  sensor  and  scatterers  are  located  within  in  a  homogeneous,  isotropic, 
locally  linear  dispersive  material  with  a  known  relative  dielectric  permittivity  e(u;). 
Under  the  assumptions  of  a  scalar  model  of  wave  propagation  and  single-scattering, 
the  electric  field  component  of  the  scattered  field  is  given  by 

r  e-iuj(t-2n(u))\rs,y\/c0) 

(!)  Esc(x,  t,s)  =  /  — — — - -2 — u2P(uj)dujf(y)dy, 

Jc  (167r3|rSiy|) 

where  P(uj)  is  the  spectrum  of  the  transmitted  pulse,  T  is  a  modified  target  that 
accounts  for  both  the  target  reflectivity  and  non-planar  surface  area,  y  =  (2/1, 2/2)  is 
a  two  dimensional  vector  position,  and  C  is  a  Bromwich  contour  in  the  upper  half 
plane.  The  complex-valued  index  of  refraction  of  the  dispersive  material  is  n( u>)  = 
a/ e(cu),  the  flight  path  q(s)  is  parameterized  by  the  slow  time  s,  and  rs  y  =  |\&(y)  — 
q(s)  |  is  the  distance  between  the  antenna  and  the  target.  Here,  it  is  assumed  that 
the  targets  are  stationary  and  are  comprised  of  linear  materials,  the  antenna  is  an 
isotropic  point  source,  the  same  antenna  is  used  to  transmit  and  receive,  the  antenna’s 
position  remains  fixed  from  transmit  to  receive,  scattering  occurs  at  a  known  surface 
^(x1,x2)  =  [xi,Xi,ip(xi,X2)]T ,  and  the  target  dispersion  in  known. 

It  is  assumed  that  the  received  signal  is  corrupted  by  additive  white  noise  r/(s,t). 
The  image  is  formed  by  applying  a  filter  Q  to  the  noisy  data  and  back  projecting 
using  only  the  real  part  of  the  refractive  index  nr(u>)  in  order  to  account  for  phase 
delay 

/( z)  =  J  /c°)Q(c/,s,z)du/ 

(2)  •  [Esc(s,  t)  +  r/(s,  t)]  ds  dt. 
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The  dielectric  permittivity  of  the  dispersive  material  is  modeled  by  the  Fung-Ulaby 
model  [7]  for  vegetation  with  the  same  material  parameters  as  used  in  [5,  6].  The 
Fung-Ulaby  model  gives  the  relative  dielectric  permittivity  as  a  mixing  formula  for 
water  and  leaf  as 


(3)  e(u>)  =  vi€i(u )  +  (1  -  v{), 
where  ^  =  0.10  is  the  fractional  volume  of  the  leaves  and 

(4)  €i(u)  =  5.5  +  ,5'5- 

1  —  IUJT 

Here,  em  =  5  +  5 1.5  thy,.  vw  is  the  fractional  volume  of  water,  and  r  =  8ns  is  the 
relaxation  time. 


The  Scattering  Precursor.  For  a  target  T  that  is  an  isotropic  point  scatter  in¬ 
dependent  of  frequency,  each  component  of  the  scattered  electric  field  is  given  by 
Eq.  (1),  viz. 


f  p— *2u;(t— n(uj)|rSjy  |/co) 

(5)  Esc(s,t)  =  /  - — ^ - uo2P(uj)dujf(y)dy 


167t3|] 


«,yl 


(6) 


oc  J  exp 


2 1  r 


s,yi 


c0 


</>(u;,6») 


(jj2P(u))d(jj, 


where  the  complex  phase  function  (p  is  defined  as 


(7) 


(u,  0)  =  iu)  [n  (u)  —  6]  =  iuo  e(u) 


and  9  =  cot/2  |rs  y|  is  a  dimensionless  space-time  parameter. 

For  values  of  9  <  \/4.5 Vi  +  1,  the  contour  may  be  enclosed  in  the  upper  half  plane, 
which  gives  Esc(s,t )  =  0  for  t  <  (2|rSjy|/c)\/4.5u/  +  1.  For  values  of  9  >  \/4.5 v\  +  1, 
an  asymptotic  approximation  to  Esc(s,t),  valid  for  large  propagation  distances  |rs,y|, 
may  be  obtained  by  use  of  a  uniform  asymptotic  expansion  of  the  integral  representa¬ 
tion  Eq.  (1).  There  is  one  accessible  saddle  point  tvsp(9)  of  the  complex  phase  function 
<p(ui,9)  that  moves  down  along  the  positive  imaginary  axis  for  \/4.5 Vi  +  1  <  9  <  90, 
crosses  the  origin  at  the  space-time  point  9  =  90,  and  continues  down  the  negative 
imaginary  axis  approaching  the  branch  point  at  —i/r  for  9  >  90.  The  saddle  point 
usp(9)  coincides  with  the  amplitude  critical  point  uj2  appearing  in  Eq.  (1)  when  9  =  9q. 
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A  uniform  asymptotic  method  is  obtained  following  Bleistein  [8]  with  the  result 

-3/2 


Esc(s,t )  ~  a0  (  2  rs’y 
Co 


W2  71 


l«,y 


Co 


(8)  +  «i(^£ 

V  c0 

as  |rs>y  |  — >■  (X).  Here, 

/x  /cn2\  .  duj 

(9)  «„=^)pm- 

(10) 

(11) 

(12) 
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7 

Q-o  — 

(t, 

) 

Wre(0  =  \/27r  e^2/2  for  n  =  0, 1,2, - 


An  example  of  a  scattering  precursor  waveform  is  shown  in  Fig.  (1)  as  a  function  of 

9. 


Comparison  to  Previous  Work.  The  optimal  waveform  of  Varslot  et  al.  [5]  is  ob¬ 
tained  by  minimizing  the  mean-square  error  between  the  reconstructed  image  Eq.  (2) 
and  the  ideal  image 

(13)  Injz)  =  [  e'C-^UyJdJ, 

Jqz 

where  fiz  defines  the  data  collection  manifold,  the  set  of  Fourier  components  of  T 
that  are  present  in  the  data,  and  (cu,  s)  — »■  £  is  a  Stolt  change  of  variables.  The 
minimization  is  subject  to  the  constraint  that  the  total  transmitted  energy  along  the 
flight  path  7(s)  is  fixed.  The  minimization  leads  to  a  fixed  point  algorithm  that  is 
solved  numerically  to  give  the  magnitude  of  the  spectrum  of  the  optimal  waveform. 
The  minimum  phase  waveform  of  this  spectrum  then  gives  the  optimal  waveform. 
The  optimal  waveform  is  dependent  upon  the  signal-to- noise  ratio  (SNR).  Figure  (2) 
shows  the  optimal  waveforms  derived  by  Varslot,  et  al.  for  various  signal-to-noise 
ratios. 

Comparison  of  the  scattering  precursor  given  in  Fig.  (1)  with  optimal  waveforms 
for  high  SNR  in  Fig.  (2)  shows  that  the  two  waveforms  are  similar  although  they 
are  derived  in  very  different  ways.  The  scattering  precursor  is  given  by  an  analytic 
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Figure  1 .  A  scattering  precursor  waveform  for  the  Fung-Ulaby  model 
of  dielectric  permittivity. 
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Fig,  6.5.  Optimal  transmit  waveform,  as  determined  according  to  (4,8),  for  varying  levels  of 
SNR:  ^0,  20,  10,  0,  -10,  and  -20,  respectively . 


Figure  2.  Optimal  waveforms  for  various  signal  to  noise  ratios  given 
in  T.  Varslot,  J.  H.  Morales,  M.  Cheney,  Waveform  design  for  synthetic- 
aperture  radar  imaqinq  throuqh  dispersive  media,  SIAM  J.  Appl.  Math., 
71  (2011),  pp.  1780  -  1800. 
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expression  based  upon  the  material  properties  through  which  the  wave  travels  whereas 
the  optimal  waveform  requires  the  numerical  solution  of  a  minimization  problem. 


Future  Work.  In  order  to  finish  our  analysis  of  waveform  design  for  synthetic  aper¬ 
ture  radar  imaging  through  dispersive  material,  resolution  analysis  for  both  waveforms 
and  a  comparison  of  the  image  reconstructions  for  the  two  waveforms  is  needed.  This 
work  is  forthcoming. 
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